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rule that relates the nonperturbative shape functions for these weights to the shape 
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implicit in NLL resummation, and discuss its limitations. In the absence of data 
analysis for the new event shapes, we compare these predictions to the output of the 
event generator PYTHIA. 
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1. Introduction 

The hadronic final states of hard-scattering reactions encode information about the 
full range of QCD dynamics, from the short-distance production of quarks and gluons 
reflected in jet production, to the long-time process of hadronization reflected in the 
exclusive distributions of particles in the final state. Event shapes Jj], ||, which 
generalize jet cross sections, are functions of final-state momenta that measure the 
flow of energy ||. We refer to these functions as 'weights' below. Event shape weights 
are tools with which we may analyze the long distance behavior of QCD through 
jet events. Qualitatively, we expect that cross sections for states with narrow, low- 
mass jets are sensitive to long-distance dynamics, while more inclusive cross sections, 
dominated by higher-mass jets, are predicted more accurately by perturbation theory. 
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By varying the values of event shape weights we may select narrow or wide jet events, 
and thus vary the relative importance of long- and short-distance dynamics. 

Infrared safe event shapes are those that are well-defined order by order in pertur- 
bation theory. For the mean values of such event shapes, or for values of the weight 
where higher order perturbative corrections are relatively small, nonperturbative 
effects enter as additive power corrections. The analysis of perturbation theory sug- 
gests that these are typically integer powers ofl/Q, with Q the overall center of mass 
(cm.) energy in annihilation processes or the momentum transfer in deep-inelastic 
scattering || [7|, §, |9|, [H| [11], [12| . The universality properties of power corrections, 
also inferred from perturbation theory, can be used to provide measurements of the 
strong coupling [13|, [14| . As we approach the narrow jet limit, however, perturbative 
cross sections generically develop large logarithmic corrections. 

When the weights of event shapes are symmetric in the phase space of final 
state particles, their large logarithmic corrections may be resummed in the limit of 
narrow jets, to next-to-leading logarithm (NLL) JT5], [R| [T7|] and beyond jrS[] in e + e~ 
annihilation. The analogous but even more challenging analysis has been carried out 
for jet events in deep-inelastic scattering |1| ^0|. For many such quantities in e + e~ 
annihilation, it is possible to factorize perturbative and nonperturbative corrections 
in a well-defined fashion and to study the latter phenomenologically. 

Familiar examples include the thrust, T [|IJ, jet broadening, B and the heavy 
jet mass, pu 0- Short distance effects are treated in resummed perturbation theory 
and control the leading-power behavior in Q including logarithms [[1|, [16], [T7|, |T8[ |. 
Long-distance effects still enter as power corrections, whose form may be inferred 
from perturbation theory in a manner analogous to power corrections to the mean 
values of event shapes f| ||. An example is 1 — T, with T the thrust in e + e~ 
annihilation, for which the first power correction near T = 1 is 1/[(1 — T)Q\. This 
correction grows rapidly as T — > 1, but it exponentiates ||, |{|. The effect is to 
shift the perturbative distribution B 0, which allows additional determinations of 



the strong coupling in conjunction with nonperturbative parameters [14]. For small 



1 — T, not only the first power correction, but all corrections of the form 1/[(1— T)Q] n 
may be organized into shape functions [21], |22j in a manner that we will review below. 
Considerable progress has been made in determining shape functions for the 



thrust, C-parameter and heavy jet mass |23|, ^2] |25], [26|, fT7|| . It is clearly desirable, 
however, to widen the class of event shapes to allow systematic studies of the ap- 
proach to nonperturbative dynamics in infrared safe observables. 

In this paper, we attempt to pursue such an approach, using the class of event 



shapes introduced in ||18| , which includes the thrust and jet broadening as special 
cases. As in [IS], we will restrict ourselves to event shapes that vanish in the two-jet 
limit for e + e~ annihilation, but extensions to quantities that probe multi-jet structure 
are possible |2"Bfl , and we will discuss the differences briefly below. 

We will follow Refs. || [^, |2T| in deducing the form of power corrections from the 
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resummed event shape weights at next-to-leading logarithm |1| in terms of shape 
functions. For the class of weights that we consider, we will identify an approximate 
scaling relation for their nonperturbative shape functions. We argue that this scaling 
reflects the boost invariance of the underlying soft dynamics. We will work for the 
most part directly in a transform space, where the resummation of large perturbative 
as well as power corrections is most transparent. 

Because we know of no experimental data analysis for this class of observables, 
we compare the predictions of our scaling rule to the event generator P YTHIA [59] . 
As we shall see, PYTHIA reproduces the scaling over a range of parameters. We 
will briefly discuss the dynamical effects that are neglected in the derivation of our 
scaling law, and what we may conclude from the event generator results. 

The following section reintroduces the event shapes of Ref. ||18|| , and in Section 
3 we review their resummation, including their numerical evaluation in perturbation 
theory, following the methods of Refs. fli~5 , |l~6[1 . We derive the scaling rule and 
describe its interpretation in terms of boost invariance in Section 4. In Section 5 we 
discuss its comparison with PYTHIA [^9| for shifts in the distribution peaks and for 
moments of the shape functions. We close with a summary. A brief appendix gives 
some details on the evaluation of integrals for the resummed distributions. 



2. The Class of Event Shapes 

We consider final states N characterized by two nearly back-to-back jets, 

e+ + e- -> J 1 (N) + J 2 (N), (2.1) 

at center of mass energy Q ^> A QCD . For each such final state, a thrust axis is 
determined as the unit vector n that maximizes the the thrust T ]]]: 

E \Pi ■ n\ 

T(A)=max alli g! — t . (2.2) 

* E \Pj\ 

all jeN 



The class of event shapes r a of Ref. |18| is defined in terms of this axis. Parameter a 
is adjustable, — oo < a < 2, and allows us to study various event shapes within the 
same formalism. It helps to control the approach to the two-jet limit. 
The weight functions defined in Ref. [pl| for a state N are 

r a {N) = i J2 ( sin ^) a - I cos ^l) 1_a > ( 2 - 3 ) 

^ all ieN 

in terms of particle energies iOi and the angles 0, of the particle momenta to the 
thrust axis. For the discussion below, we will find it more convenient to introduce a 
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form directly in terms of particle momenta 

1 

Q 



T *W = ^ £ ftTe- WM , (2.4) 



all ieJV 

where pi t is the transverse momentum relative to the thrust axis, and rji is the 
corresponding pseudorapidity, 

rji = lncot(0i/2) . (2.5) 

The case a = in Eq. ( |2.4j ) is essentially 1— T, with T the thrust while a = 1 is the jet 
broadening. A similar weight function with a non-integer power has been discussed 
in a related context for 2 > a > 1 in M . The two forms for the weight functions, Eqs. 
( |2.3|) and fl2.4|) are equivalent for massless particles. We will use Q2.4|) , in terms of 
momenta and pseudorapidities, in the analysis of nonperturbative corrections below. 
We will recall later the difference between the two forms of the event shapes when 



applied to massive hadrons 30 



Any choice a < 2 in ( |2.4|) specifies an infrared safe event shape variable, because 
the contribution of any particle i to the event shape behaves as 6i 2 ~ a in the collinear 
limit, 9i = cos -1 \hi -h\ — > 0. However, the resummed formula given below from Ref. 
[ [iljfl is valid only for a < 1. For a > 1 recoil effects have to be taken into account, at 
least beyond the level of leading logarithm [O]. As a — > 2 the weight vanishes only 
very slowly for Bi — > 0, and at fixed r tt , the jets become very narrow. On the other 
hand, the limit a — ► — oo corresponds to the total cross section. 

The differential cross section for such dijet events at fixed values of r a is given 

by 

d ~^f 1 = \M(N)\ 2 S(r a -r a (N)), (2.6) 

where we sum over all final states N that contribute to the weighted event, and 
where M(N) denotes the corresponding amplitude for e + e~ — > N. 

Since we are investigating two-jet cross sections, we fix the constant r a to be 
much less than unity: 

< r a < 1 . (2.7) 
For small r a , the cross section ( |2.6| ) has corrections in ln(l/r a ), which have been 



organized in Ref. [18]. In the following we will quote the result of the resummation 
of large logarithms of r a in Laplace moment space. The Laplace transform of the 
cross section fl2.6| ) is given by 

Logarithms of l/r a are transformed to logarithms of v. For large u, dependence 
on the upper limit in the r a integral is exponentially suppressed. Here and below 
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quantities with tildes are the transforms in r a , and quantities without tildes denote 
untransformed functions. Our results below are valid in the region where In v is much 



larger than \a\ |L8 



3. The Resummed Cross Section at NLL 

3.1 The resummed cross section in moment space 

The NLL resummed cross section Q2.8Q for a < 1 in moment space can be written as 



IS 



Ctot 



a (u, Q, a) = exp < 2 



du 
u 



uQ 2 



dp 



T 



Pt 



A(a s (p T )) (e-" 1 " ^) - l) 



u 2 Q 2 

+\B (a s (Vu-Q)) ( e -«^) 2/(2 - a) - 1 



(3.1) 



where Jiv^ Q, a) is a factorized function associated with each jet. The resummation 
is in terms of anomalous dimensions A(a s ) and B(a s ), which have finite expansions 
in the running coupling, 



(3.2) 



n=l 



and similarly for B. To NLL they are specified by the well-known coefficients, 



a« = c F , 

A(2) = \ cf 



1 18 6 J 9 ; 



- — Cp . 

2 ' 



(3.3) 
(3.4) 

(3.5) 



independent of a. Cp and Ca are the Casimir charges of the fundamental and adjoint 
representation of SU(iV c ), respectively, Nf denotes the number of flavors, and T F = 
1/2 is the usual normalization of the generators of the fundamental representation. 
Eq. (|3.1|) reproduces the NLL resummed thrust cross section [15], |16| when a = 0. 



3.2 Inversion of the transform 

As it stands, the resummed cross section ( |3.1|) is ambiguous, because of the singularity 
of the perturbative running coupling. To define the integrals in ( |3.1| ) and to invert 
the transformed cross section from moment space back to r a , we will follow the 
method of Ref. ||16|| . In this approach, we avoid the singularities of the perturbative 
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running coupling by reexpressing the running coupling in terms of the coupling at a 
hard scale, and by performing the resulting integrals in the exponents of Eq. (|3.1| ) 
to NLL. Explicit expressions for the cross section, and hence the jet functions J7", 
are given in Appendix |A|. The inversion is then also carried out to NLL. In the 
final expressions, the singularities of the running coupling are only manifested as 
singularities at small values of r a . 

As in Ref. [|16||, we work with the integrated cross section, also called the radiator, 



R(r a , Q) 



Ctot 



dr. 



i da(r' a , Q) 



dr' a 



(3.6) 



The radiator can be found directly from the jet functions, J{y, Q, a) in transform 
space by 



R{r a , Q) 



1 

2vri 
1 

2m 



c 



due UTa R{u,Q,a) 
-e^[J(u,Q,a)} 2 



(3.7) 



c 



The contour C lies in the complex plane to the right of all singularities of the inte- 
grand, and the first equality defines the radiator in transform space. The differential 
cross section ( |2.6D is easily obtained once we have an explicit form for the radiator, 

1 da(r a ,Q) 1 d 



-R(r a ,Q). 



(3.f 



otot dr a r a dhiT a 

To perform the integral in Eq. (|3.7|) , we Taylor expand the resummed exponent 
with respect to lnv around lni/ = ln(l/r a ), because the functions g\ and g 2 in J 
(see Appendix [A]) vary more slowly with v than vr a [16]. At NLL accuracy we can 
neglect all derivatives higher than first order in the Taylor series of the exponent. 
Performing the integral is then straightforward, using 



2m J c r(l- 7 ) 



(3.9) 



In this way, we find 



R(r a ,Q) 



exp |21n (i) 9l (x, a) + 2g 2 (x, a) + 2(2 - a)x 2 In {^j g[(x, a) j 



1 - 2g 1 (x,a) - 2xg[(x,a) 



(3.10) 



The functions g iy i = 1, 2, and g[ are given in Eqs. ( |A.6|) , ( |A.7|) , and ( |A.8|) of Appendix 
^ in terms of the variable 



X = rrrz r m 



7T 2 (2 -a) \t ( 



(3.11) 
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with (3q the first coefficient of the QCD beta function. The explicit formulas in 
the appendix show that the gi(x,a) have logarithmic singularities at x = 1 and 
x = 1/(2 — a), which are the manifestations of the singularities of the perturbative 
running coupling for this NLL evaluation. Although well defined for most values of 
r a , the resummed cross section remains ambiguous beyond NLL. 



3.3 Matching and numerical evaluation 

A realistic evaluation of the resummed cross section requires matching to fixed-order 
calculations. While the resummed predictions are reliable for small values of r OJ fixed- 
order contributions are more accurate at larger r a . For matching to NLL accuracy, 
it is necessary to know the fixed order contributions up to 0(a 2 ). These can be 
calculated numerically using, for example, the program EVENT2 f32]| . 

In the following we will use log R matching JjlJ. Other matching schemes are 
possible, and differ from this formula at order ot z s and NNLL. The radiator at NLL, 
evaluated at the scale \i = Q, is computed as 



lni? NLL (r a , Q) = lni? rcsum (r a , Q) + lni? fl 



,Q)-lnR:Z m (r a ,Q) • (3.12) 



Here the first term on the right, lni? resum , is the logarithm of the resummed radiator, 
which can be read off of Eq. (|3.10|) . In -R flxed is the logarithm of the fixed order radiator 
calculated with EVENT2 at NLO, expanded to order a 2 , 



fixed fixed + 



7T 



<*s(Q) 



7T 



R 



(2) 



d(1) 



(3.13) 



where the R^ ed , i = 1,2 are the first- and second-order parts, respectively, of the 
radiator given by EVENT2. The last term on the right of Eq. ( ^TT2|) , In i?^ s p um (r Q , Q), 
is the logarithm of the resummed radiator expanded to order a 2 , which needs to be 
subtracted in order to avoid double counting. This contribution is found from Eq. 
( ETTOD with Eqs. flOfi , (|Ap , and JOj) of Appendix 0: 



In i? exp 

rcsum 



Ta: Q) 



+ 



a s (Q) 

71 

a s (Q) 

71 



C7 L1 ((,)ln( lj +G 12 (a)ln 2 \J- 



2 r 



(V 2 o(r/)hr ( y\ +G 23 (a)lnMl 



(3.14) 



the functions Gij are listed in the appendix. Finally, the physical requirement that 
the cross section vanishes beyond the upper kinematic boundary, t™, imposes the 
following constraints on the matched expression (|3.12|) , 



R 



NLL /_max 



Q) 



dR™ (r a ,Q) 
dr n . 



0. 



(3.15) 
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This is achieved by the replacement 



16 



1 



1 



1 



(3.16) 



in the resummed terms, since the fixed order coefficients satisfy Eqs. ( p. 151 ) order by 
order. r™ ax depends on the number of final-state particles, corresponding to the order 
in a s up to which the fixed order contributions are evaluated. The replacement Q3.16] ) 
suppresses terms at large r a higher than this order that are generated inaccurately 
by the resummed contribution. r™ ax at the leading logarithmic level is given by the 
limiting configuration of three final-state particles [[35], 



T 



,LL 



V3 a 



(3.17) 



The limit from four particles in the final state which gives the upper kinematic 
boundary at NLL can easily be determined from the output of EVENT2. 

Eq. ( |3.10| ) for the radiator and hence Q3.8| ) for the cross sections, is applicable, 
with power and subleading logarithmic corrections for values of r a away from the 
end-point region, where the variable x in Eq. ( |3.11| ) becomes of order unity. That is, 
we require [/3q/ (2n)]a s ln(l/r a ) < 1, or equivalently r a > A QCD /Q. For r a ~ A 
non-perturbative corrections become dominant. In this range of r a , we "freeze 
perturbative contribution to the cross section, following [2 



< >cd/Q 
the 



23 



Rpt {r a , Q,k) = 9 



K 

t„ e 

Q 



-IE 



i T a, Q) 



-~1E 



- Q 



(3.18) 

T2|), and where k is a nonperturba- 



where i? NLL is evaluated according to (|3.7| ) and 
tive cutoff. Dynamics below the scale k will be incorporated into the nonperturbative 
corrections, in a manner we will discuss below. For our numerical studies we pick 



0.75 GeV, as in 1221. 



4. The Scaling Rule 

We are now ready to derive the scaling relation for nonperturbative shape functions. 
We show first how the rule is implied by the resummed cross section that we have 
just described, and go on to interpret the physical content of the scaling. 

4.1 From resummations to shape functions 



Following Ref. [21|, we identify the power structure of nonperturbative corrections by 
a direct expansion of the integrand in the resummed exponent at momentum scales 
below an infrared factorization scale, k. Although this scale need not be exactly 
the same as the scale in Eq. (|3.18|) at which the radiator is frozen, they are closely 
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related, and we will use the same symbol for both. We thus rewrite Eq. (|3.1|) as the 
sum of a perturbative term, summarizing all pt > k, and a soft term, containing the 
nonperturbative physics of strong coupling. This corresponds to Pt < K in the term 
with anomalous dimension A, and uQ 2 < k 2 in the term with B. Exchanging the 
order of integration for the A term, we find 

lni?(z/, Q, a) = In Rpt{v, Q, k, a) 

i 

Vt 



- p T /Q 



+2 f^fA(a s (p T )) [ ^ ( e -^(Pr/Qr _ x 
J Pt J u \ 



p 2 T /Q 2 

* 2 /Q 2 



+ J (a.cvaj)) (e-^ w -> - l) 



In Rpt(v, Q, k, a) 



2 

,, x ri(l-a) 
1-1? 



+T 4,£i(-i)7f^.(-))i-(t 

n — 1 n L 





v2n/(2-a) r2/Q2 



+ E( _ ir W!)_Jl / ^ B(as( ^ Q)) 

n=l „ 



= ln J R PT (z/,Q, K,a) + ln/ ajNP 

(41) 

In the second equality we have expanded the exponentials and integrated over u. 
In the last equality we have introduced the logarithm of the shape function, as an 
expansion in powers of v/Q, 



ln Mi- K ) = T^5^(-C> 

- ^ E w«> (-1) 



2 

T 



Pt 



(4.2) 



We neglect the terms (pt/Q)^ 1 ^ i n Eq. Q4.1|) for the expansion with A, and the 
entire B expansion, as indicated. These terms are suppressed by additional (frac- 
tional) powers of Q only, of course, for a < 1. This is the same restriction to the 



NLL resummation formula, Eq. ( |3.1| ) that follows from considerations of recoil 



Given this approximation, we find the simple result that the only dependence on a 
is through an overall factor 1/(1 — a). 
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Equation (|4.2|) immediately leads to the following scaling relation between shape 
functions for different values of a: 

ln/ - NP (i' K ) = r^ ln/6 ' NP (i' K ) • (4 - 3) 

We can now use the shape function determined by the thrust (a = 0) to predict 
shape functions for any value of a < 1. Of course, for a > 0, the relative suppression 
of the neglected terms in Eq. fl4.1| ) is by fractional powers, and we will apply our 
results only to a < below. It is worth noting that for 2 > a > 1, where the event 
shape weight vanishes very slowly in the collinear limit, the fractional powers in ( |4.1| ) 
dominate the integer powers, as observed in Ref. 0. This, however, is outside the 
range of a to which our formalism applies at NLL [IF 



4.2 Interpretation of the scaling 

The assumptions that entered our result (fOl ) are relatively standard: that nonper- 
turbative parameters may be identified with moments of the running coupling with 
respect to its scale, or more generally moments of the anomalous dimension A(a s ). 
Our analysis does not predict the values of the nonperturbative parameters A n aside 
from their overall a-dependence, and is consistent with the absence of even power 
corrections, as argued in Refs. PSL I2J1 E7|. It is also based on the specific form of 



NLL resummation. Despite these limitations, we can identify the physical origin of 
the scaling. We have seen that the leading powers come entirely from the A(a s ) 
term in the exponent of the resummed NLL cross section. This contribution can be 
derived directly from the eikonal approximation, in which all soft gluons are emitted 
from light-like Wilson lines. 

To be specific, we can define the eikonal cross section by 



r a (N)) 6{Q-E{N)) , (4.4) 



N 



where r a (N) is the weight computed according to Eq. (pO] ) for state N, and E(N) 
is the total energy of the particles in N. The operator U(0) is defined in terms 



of Wilson lines, $p[oo, y] = Pexp ig r o °° da fit 1 A}!' (ap + y)\, where (/) labels the 
representation (quark or antiquark in our case) and where P indicates color ordering 
along the corresponding light-like direction (3. In this notation, we define 

U(0)=T [4f[oo,0]<?f[oo,0]] . (4.5) 

with f3f — /9| = and /3i ■ {3% — 1, and with T the time ordering operator. We work 
in a frame where j3\ and {3% are back-to-back. 

Diagrammatic rules for constructing the perturbative expansion of Eq. ( f4.4| ) may 
be found, for example, in Ref. and the exponentiation properties applicable to 



- 10 - 



logarithmic corrections for eikonal event shapes in Refs. |34], |35]. Here we only need 
to stress that matrix elements computed from the product of back-to-back Wilson 
lines in Eq. ( [4.4]) are boost invariant along the axis defined by f3\ and {3% and are also 
invariant under the scaling of all final state momenta (pi — > Xpi) up to renormal- 
ization effects. Frame dependence enters the eikonal event shape distributions only 
through fixing the shapes themselves, which distinguish between positive and nega- 
tive rapidities (see Eq. (|2.4j)), and through ratios of the cutoff in the overall energy, 
Q, to the renormalization scale. In particular, when the decay products of a particle 
are emitted into both hemispheres [p4| , |36| , the event shape receives a noninvariant 
contribution, even if the amplitude to produce the virtual particle is boost invariant. 

These features of the eikonal cross section are illustrated by the exponent of 
the resummed perturbative cross section in Eq. (|3.1|), which reflects single-particle 



kinematics in the exponent. Comparing the definition of the weight in Eq. (2^) with 
the NLL resummation (|3.1|) , we can change variables from u to pseudorapidity (|2.5| ) , 
by it = (pr/Q) ex p( — M), du/u = — (\r}\/rj) drj. The boost and scale invariance of the 
single-particle emission cross section in eikonal approximation then become manifest. 

In general, the interplay between the limits of integration and the transform 
links the rapidity and transverse momentum integrals in Eq. (|3.1|) , to produce the 
nontrivial a-dependence of the functions gi in Eq. ( |3.10| ), given explicitly in the 
Appendix. When the transverse momentum is limited by the cutoff k of Eq. ( [4.1| ), 
however, the dependence on the remaining limits of integration is simplified. To 
see this, we expand the exponential of the A term in the first equality of ( |4.1| ). For 
Pt < k, with k < Q, we have l/[n(l — a)] < ln(Q/pr) for all n > 1, and the rj integral 
of the nth term in the expansion is cut off by the exponential exp[— n(l — a)\r]\}. 
Approximating this factor by unity when the exponent is smaller than one, and by 
zero elsewhere, each integral becomes a simple measure of the rapidity range over 
which the single-particle weight is negligible, 

k 2 HQ/pt] 

^fA (a s (p T )) J drj (e-(«PT/Q>-^^ _ i 
i) ^ T o 



oo 1 

E - 



n=l 



k z [n(l-a)]- 1 







By comparing the right-hand side here to the the right-hand side of the second equal- 
ity of Eq. (|4.1| ), we see that corrections to ( |4.6j ) are precisely the power-suppressed 
terms that we have neglected for a < 1. The overall decrease with 1/(1 — a) thus 
reflects the shrinking rapidity range that is available in each term due to the expo- 
nential of the weight. The contribution to each power correction from the region 
of strong coupling is boost invariant, but the available rapidity range decreases uni- 
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formly in (1 — a) -1 . The observed scaling, then, is a reflection of the underlying 
boost invariance of the nonperturbative dynamics. 

Alternatively, we may consider the scaling with (1 — a) -1 as a test of the use of 
NLL resummation to determine shape functions, and of the extent to which boost 
invariant dynamics dominates the differential distributions. We will return to this 
viewpoint at the end of the following section, after testing the scaling rule with the 
event generator PYTHIA. 



5. Tests of Scaling 

At present, we know of no data analysis that would provide experimental tests of our 
ability to predict nonperturbative contributions for the class of event shapes we are 
considering, aside from the thrust. In the absence of such an analysis, we will rely 



on the event generator PYTHIA |29j as a stand-in for experiment. In this section we 
will describe two tests of our scaling rule. These are, of course, tests of consistency 
with PYTHIA, not with nature. First we discuss our predictions for the shifts || [| 
in the resummed perturbative distributions due to the first power correction {yjQ\ 
and then the full scaling associated with the summation of all powers in vjQ. 

For convenience, we evaluate the weight functions using Eq. (|2.4|) , in terms of 
particle momenta, using the output of PYTHIA for the thrust axis, which is also 
computed in terms of particle three-momenta. The derivation of the scaling rule, Eqs. 
( fll|) - (fO|) does not distinguish between the weight functions expressed in terms 
of energy ( |2.3| ) and in terms of momentum ( |2.4| ), although with massive hadrons in 
the final state, the values of r a can be different. We therefore expect the scaling to 
hold only if the weight functions are computed consistently in terms of energies or 
momenta. 

5.1 Shifts of the distributions 

Let us first study the shifts of the distributions. Retaining only the term with 
n = 1 in ( |4.1|) , and suppressing the dependence on k, we obtain the cross section in 
momentum space from Eq. fl3.7|) , 



Rfa, Q) = 7T— [ -eK^-i^) [J„( V , Q, a)] 2 . (5.1) 

In this approximation the integrated cross section is shifted || || to the right from 
the perturbatively calculated spectrum by an amount 

Ar(a,Q) = -^—^. (5.2) 
1 — a Q 



To the same approximation, this also holds for the differential cross section ( |3.8| ) for 
r a not too small. 
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Here we study the shifts of the peaks, At p p5| . From Eq. 
shifts of the peaks for different values of a multiplied by (1 — 
measured at the same scale Q: 



we infer that the 
are the same when 



(1 - a) Ar p (a,Q) 



b) Ar p (b,Q). 



(5.3) 



Let us compare this prediction, valid at the partonic level, to the corresponding 
cross sections computed by PYTHIA |p9[j , version 6.215 |37]], at the hadronic level, 



using PYTHIA's implementation of the string fragmentation model ||38|| . We use the 
default settings of the program. The string picture seems to model the hadronization 
process successfully, as a comparison with recent data for the thrust shows (see e. g. 
[Hi)- As noted above, for comparison with PYTHIA we use the definition of our 
event shapes in terms of three- momenta, ( |2.4|) , for both the partonic and the hadronic 
level. Other prescriptions are of course possible |3(| . 

We see from Fig. [l] that the shifts computed this way obey the scaling quite well 
in the range between a = and a = —1. For larger values of \a\, the peaks of the 
cross sections move into regions where r a ~ A QCU /Q, and where our analysis is not 
reliable. 



0.020 



0.015- 



Q. 

co 0.010 



0.005 



0.000 



-1.0 



-0.8 



-0.6 



-0.4 



-0.2 



0.0 



Figure 1: Shifts of the peaks Ar p (a,Q = 91 GeV) of the distributions (l/cr tot ) da/dr a 
between NLL partonic resummed predictions and hadronic cross sections computed using 
PYTHIA with string fragmentation. The result is multiplied by (1 — a). The error bars 
are estimated from the uncertainty of the NLL resummed calculation and the output of 



PYTHIA. The shaded band is the shift of the peak for the thrust determined in [14] between 
resummed predictions and experimental data. 
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5.2 Shape functions 



Now let us turn to the full shape functions. We will test the scaling directly in 
moment space, where a parameterization in terms of the coefficients A n of Eq. (|4.2|) 



is most straightforward. Using Eq. ( |4.1| ), we have 



/a,NP ( q, « 



V 



R(u, Q, a) 
Rpt{v,Q,k,(i) 



+ o 



Q 



l-a 



(5.4) 



This gives moment space expressions for the shape functions directly, given exper- 
imental or other input for R(v, Q, a) in the numerator and resummed perturbation 
theory (in our case to NLL) for Rpt(v, Q, k, a) in the denominator. As above, we use 
the output of PYTHIA for the numerator. For the denominator, we use the method 



described above in Sec. 3.3 



The results are shown in Fig. || at Q = 91 GeV and in Fig. ^| at Q = 35 GeV. 
We begin by computing the shape functions for a = 0, —0.25, —0.5 at Q — 91 GeV 
and Q = 35 GeV directly from Eq. ( |5.4j ). These are the solid curves in the figures. 
The dotted curves show the predictions found by simply scaling the a = curve in 
each case according to Eq. ( |4.3| ), or equivalently 



/a,NP ( q, K 



V 



/o,NP ( q, K 



V 



(5.5) 



The perturbative radiator defined in Eq. (|3.18|) , and thus the ratio in ( |5.5| ), as plotted 
in Figs. ^ and |3|, are fairly dependent on the value of the cutoff k. This dependence 
compensates of course between perturbative and nonperturbative contributions to 
the full radiator (|0|). 

As in the case of the shifts, we see that within the range of parameter a consid- 
ered, the scaling works well. It is worth pointing out that a similar scaling fails by a 
relatively large amount for the perturbative cross section itself. We have restricted 
ourselves to a minimum value of a = —0.5 to keep an acceptable numerical accuracy. 

We note that the logarithms of the curves in both figures depend fairly linearly on 
v for relatively small u, indicating that in this range the linear, vjQ term dominates. 
This is consistent with the result above for the shifts, which is based only on the 
vjQ correction. To fit the curves at larger u, however, higher powers in vjQ are 
necessary. While we have not attempted a fit of the a = case (thrust) at different 
energies our reasoning is consistent with any determination of the coefficients A n (/t) 

in Eq. (O) pa papa]. 



We emphasize that the agreement of our scaling rule with PYTHIA may only 
mean that the output of PYTHIA shares some of the properties that go into the 
derivation of the rule. 



- 14 - 



1.0 





Figure 2: The ratios, Eq. (5.4), between the output of PYTHIA and the NLL resummed 
predictions at a = 0,-1/4,-1/2 at cm. energy Q = 91 GeV. The solid lines are the 
directly computed ratios, the dotted lines is the scaled a = curve, according to Eq. (|5.5|). 



5.3 Scaling violations 



wc 



In the light of the relation between the scaling rule and boost invariance, Sec. \L2 
can understand why PYTHIA respects the rule over a moderate range of parameter 
a. It is natural to think of jet fragmentation as dynamically boost-invariant along the 
jet axis, while boost invariance need not be respected for coherent interjet radiation, 
which depends on the relative directions of the jets. For moderate values of a, the 
shape functions may follow the scaling in our numerical tests because in PYTHIA the 
corresponding event weights are dominated by particles created from boost-invariant 
dynamics, and for which correlations between the jet hemispheres (for example, due 
to decays into opposite hemispheres) are negligible. 

Beyond NLL, however, the general resummation for the two-jet event shapes 



of Eq. (|2.4j) involves coherent interjet radiation | l8fl , which we expect to produce 
correlations between the jet hemispheres. The neglect of such correlations is related 
to the "inclusive" approximation discussed in Ref. in which the effect on the 
weight function from off-shell gluons that split into particles that move into different 
hemispheres is suppressed. 

More generally, suppose that the correlation between energy flow is enhanced 
relative to NLL perturbation theory for "short-range" rapidity intervals, less than 
some constant At]o- Then, because of correlations between radiation in different 



hemispheres, the scaling implied by Eq. (O) can hold only when the rapidity range 
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1.0 





Figure 3: The ratios, Eq. (5.4), between the output of PYTHIA and the NLL resummed 
predictions at a = 0, —1/4, —1/2 at cm. energy Q = 35 GeV, solid and dotted lines as in 
Fig- | 



for each term on the right-hand side is much larger than the constant, [n(l — a)] ^> 
Ai]q. For any Arj , this condition is violated for n large enough. Thus, as (1 — a) 
increases, the scaling is violated, first by high powers of v/Q, and eventually by 
lower powers. Only to the extent that enhanced short-range correlations are fully 
negligible can we expect that the scaling holds for the full shape function. Pushing 
our analysis to larger values of (I — a) should eventually uncover correlations between 
hemispheres, even in PYTHIA, perhaps associated with the string breaking picture 
of hadronization. In general, we would expect physical moments to decrease less 
rapidly than (1 — a) -1 in the presence of positive correlations between radiation in 
different hemispheres [|4[ . 



The extension of our analysis to multijet events []5S| should be relatively straight- 
forward. In this case, as in the two-jet limit of e + e~ annihilation, it is necessary to 
determine each jet axis by a thrust-like condition, which is relatively insensitive to 
recoil effects [17, 18], to insure resummation at NLL. Such cross sections still factorize 
perturbatively, but now into a function involving coherent interjet radiation as well as 
jet functions, even at NLL. The jet functions, but not the soft function, should obey 
a scaling relation like Eq. (L3). Tests of this scaling for multi-jet shape functions 



could be an indirect way to estimate the significance of coherence effects for interjet 



radiation. Extensions to deep-inelastic scattering |]T9| and hadronic scattering p8 
may also be possible. 
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6. Conclusions 



We have derived a scaling property that relates nonperturbative shape functions 



within the class of event shapes, including the thrust, introduced in fll8fl . We have 
seen that the resulting predictions match the output of the event generator PYTHIA 
over a range of the relevant parameter that defines the event shape. This analysis 
is based on the NLL resummed cross section for these event shapes, which neglects 
correlations between hemispheres. The comparison of these predictions to actual 
data should shed light on universality properties of nonperturbative corrections, and 
their relationship to determinations of a s . We hope that this comparison is still 
possible for archived data from LEP. 

We have argued that the scaling rule derived from the NLL cross section is more 
generally dependent on the boost invariance of strong coupling dynamics. In general, 
the scaling must also fail at some level due to correlations between hemispheres, 
associated for example with decays [|4|, |36| . We might expect such effects to become 
more important for a <C — 1, since the corresponding weight functions are sensitive 
primarily to radiation near the boundary between the hemispheres. In any case, we 
hope that the example studied above shows that event shapes can be designed to 
probe specific aspects of nonperturbative QCD dynamics. 
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A. Explicit Expressions for the Cross Section in Transform 
Space 



For the evaluation of the integrals in Eq. (|3.1|) we use the running coupling at renor- 
malization scale /i in terms of the coupling a evaluated at Q/2, expanded for use at 
NLL accuracy, 

'?) 



a = a* 



(a.i; 



a 



1 



a 



47rA)l + ^aln^ 



ln|1 + ^ ln |) + ... 



where the coefficients (3 and Pi are given by 



(A.2) 



(A.3) 
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34 2D 
A = y C\ - -C A T F N f - ACpTpNf. 



(A.4) 



The term with Pi is only necessary for the integral containing at NLL. 

Inserting the expansion of dA~T|) into Eq. (g^), the integrals are done to NLL 
accuracy in terms of elementary functions by using the replacement (exp[— y] — 1) — > 
—9 (y — e~ lE ), which is accurate to NLL. Here 7^ is the Euler constant. The result 
of this procedure is, 

1 

-a(v,Q,a) = [ J(u,Q,a)} 2 



Otot 



exp 1 2 ln(v) 5-1 
+ 2 
+ 2# 2 



/ A) «s 


Mm,, 


•) 


V2tt2 


— a 






In 2 ^ In 


/2/i 


2 - a 







2/i\ , / Po a s (jj) 



9i 



2vr 2 - a 



In za a 



2vr 2 - a 



In*/, a) +C«ln n - 1 z/) i, 



(A.5) 



where the functions gi and #2 and g[ that resum leading and next-to-leading loga- 
rithms, respectively, are given by 



gi{x,a) 
g 2 (x,a) 



4 1 1 
Pol — ax 



1 



2 - a 



x) ln(l - (2 



ax 



'1- x) In(l-x) 



(A.6) 



|-BW ln(l -x)- ^L-^—A^ [(2 - a) ln(l - x) - ln(l - (2 - a)x)\ 
Po Po 1 — a 

- ^rlE-^—A^ [ln(l -x)- ln(l - (2 - a)x)] 
Po 1 — a 

+ — ln2—— [(2 - a) ln(l - x) - ln(l - (2 - a)ar)] 
Po 1 — a 



Pi 1 
"po 3 l-a 



[2 ln(l - (2 - a)x) - 2(2 - a) ln(l - 2) 
+ ln 2 (l - (2 - a)x) - (2 - a) ln 2 (l - a; 







Qi (x,a) = — gi(x,a) 



(A.7) 
(A.8) 



In Eq. ( A.5 ) the scale \i should be chosen of the order of the hard scale to avoid 
further large logarithms, and we choose /1 = Q. Setting \i = Q as in |Tj| cancels the 
term proportional to In 2 in Eq. (|A.7|) , and we reproduce for a = the form of [|I~6 |. 

Finally, for sake of completeness, we list the coefficients Gij that occur in the 
expansion in terms of a s of the logarithm of the resummed radiator, Eq. ( |3.14| ), 



G11 (a) 



2 - a 



(A.9) 
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G l2 (a) 
G 2 2(a) 

G 23 (a) 



2 - a 
1 1 



(2 - a) 2 36 
3 — a 1 



C F , (A.10) 
C F [48vr 2 CV + (169 - 134 a - 6 (2 - a)n 2 ) C A 



(2 - a) 2 9 



-2(11 - 10 a)N } \ 
C F {IIC A - 2N f ) . 



For a = 0, the Gy-s again reduce to those listed in [16 



(A.ll) 
(A.12) 
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